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Abstract We present a technique for automatic determination of flare-ribbon 
separation and the energy released during the course of two-ribbon flares. We 
have used chromosphcric Ha filtcrgrams and photosphcric line-of-sight magne- 
tograms to analyse flare-ribbon separation and magnetic-field structures, respec- 
tively. Flare-ribbons were first enhanced and then extracted by the technique of 
"region growing" , i.e., a morphological operator to help resolve the flare-ribbons. 
Separation of flare-ribbons was then estimated from magnetic polarity reversal 
line using an automatic technique implemented into Interactive Data Language 
(IDL™) platform. Finally, the rate of flar e-energy release was calculated us- 
ing photosphcric magnetic-field data and the corresponding separation of the 
chromosphcric Ha flare-ribbons. This method could be applied to measure the 
motion of any feature of interest (e.g., intensity, magnetic, Doppler) from a given 
point of reference. 
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1. Introduction 



Solar flares are energetic transient events, which are produced by reconnec- 
tion of magnetic-field lines at coronal heights ( |Sturrock, 1966| |Hirayama, 1974| 
|Kopp and Pneuman, 1976[ ). When a coronal flux rope loses equilibrium and trav- 
els upwards, an extreme reconnection current sheet (RCS) is formed underneath. 
The reconnection in this RCS releases most of the magnetic energy stored in the 



magnetic- field configuration ( Forbes and Priest, 1984 Lin and Forbes, 2000). Charged^ 
particles can be effectively accelerated by electric field in the RCS ( Martens and Young, 1990[ | 
Litvinenko and Somov, 1995). Some of these energetic particles, produced dur- 



ing a solar flare, gyrate around the field lines and propagate toward the un- 
derlying footpoints, precipitating at different layers of the solar atmosphere to 
produce two bright lanes, or flare- ribbons. 
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The Yohkoh Soft X-ray Telescope (SXT) observed a two-ribbon structure for 
the first time in the hard X-ray energy range above 30 keV, suggesting that elec- 
trons are accelerated in the whole system of a flare arcade ( |Masuda, Kosugi, an d Hudson, 200 ip .| 
They analyzed the motions of two hard X-ray ribbons assuming these to be the 



footpoints of reconnected loops. |Zhou, Ji, and Huang| (|2008p studied footpoint 



motion of two large solar flares, including the X10 flare of 29 October 2003 
using the UV/EUV observations by the Transition Region and Coronal Explorer 
(TRACE) and hard X-ray (HXR) data by Reuven Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI). They used the center-of-mass method to locate 
the centroids of the flare-ribbons. 

Evolution of two-ribbon flares is morphologically characterized by separation 
of two ribbons in the chromosphere, traditionally observed in Ha. This usually 
occurs during solar eruptive phenomena (eruptive filaments, flares and CMEs) 
and is believed to be the lower atmospheric signature of magnetic reconnec- 
tion progressively occurring at higher levels subsequent to energy release in the 
corona. As the magnetic-field lines at coronal heights progressively reconnect, 
their footprints move farther out from the already reconnected field lines. This 
successive reconnection process causes the apparent increase of the separation of 
flare-ribbons. As a result of magnetic reconnection, loop arcades are produced 
below the reconnection site. These loops are initially very hot and visible in EUV 
(c/., FigurcQJi. b), and only the footpoints of the loops are observed in Ha (c/., 
Figure [T};). The two sets of footpoints at each side of the loop arcade appear 
as two elongated ribbons. As the loops cool down, the whole or parts of loops 
may be temporarily visible in Ha. These are called Ha (post-) flare loops. The 
separation velocity of flare-ribbons depends on the reconnection rate of magnetic- 
field lines, indicating a close relationship between flare-ribbon separation and 
energy release. 

Two well-known reconnection models are the Sweet-Parker and the Petschek 
models. Reconnection rate is one of the most important physical quantities in 
reconnection physics. ISweetl(|1958p and IParkerl(|1957p predicted low reconnec- 
tion rate due to high Reynold number (R m ss 10 8 - 10 12 ) in the solar corona, and 
therefore the magnetic-field lines are frozen. As a result, the Sweet-Parker model 
is inefficient in producing fast reconnection, and conversion of magnetic energy 
to plasma energy takes place at a slow rate. However in a thin current sheet 
in the solar atmosphere, R m could be very small (« 1) so the magnetic-field 
would not be frozen and could slip through plasma while the magnetic energy is 
converted to heat in short time scale QPetschek, 1964| ) . This model may explain 
the energy-release rate and time of solar flares, and it is now widely accepted 
that the main driver of the solar flare is fast coronal magnetic reconnection of 
Petschek type. 

The main requirement in understanding the flare reconnection rate and en- 
ergy release is to detect and obtain flare characterization by image processing 
and pattern recognition techniques. Several methods have been proposed for 
automatic tracking of the apparent separation motion of two-ribbon flares. 

We have developed a technique for detection of flare, its characterization and 
determination of various physical parameters required for calculation of the re- 
connection rate and energy released during the flare. Particular attention is given 
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Figure 1. Images (in negative) of NOAA 10486 taken on 28 October 2003 during the X17/4B 
flare: (a) TRACE UV 284 A at 11:14:12 UT, (b) TRACE UV 1600 A at 11:14:35 UT, (c) USO 
Ha at 11:14:45 UT showing post-flare loops, and (c) RHESSI HXR at 11:14:03 UT in the 
energy range 100 - 200 keV. 



to the calculation of these quantities over different parts of flare-ribbons and not 
on the flare as a whole. We have applied our technique to the extensively observed 
Xf 7/4B flare of 28 October 2003, and separation velocity of the flare-ribbons is 
estimated perpendicular to the magnetic neutral line. Velocities of different parts 
of flare-ribbons are then used to determine the energy release rate. In Section^ 
we briefly discuss the basic formalism of the problem and Section [3] describes 
the data used in this study. Data processing steps are discussed in Section |4j 
The method to determine flare-ribbon separation is discussed in Section [5] flare- 
ribbon expansion and energy release rate from different parts of the solar flare 
of 28 October 2003 are described in Section [6] Finally, the conclusions are given 
in Section [7] 
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2. The Basic Formalism 



Properties of magnetic reconnection in the corona have been related to observed 
signatures of solar flares, e.g., IForbes and Lin (2000]) . Their formulation requires 
measurement of photosphcric magnetic-fields and flare-ribbon separation speeds 
which can be used to derive two physical terms for magnetic-reconnection rates: 
the rate of magnetic-flux change involved in magnetic-reconnection in the low 
corona and the electric field inside the RCS that is generated during magnetic- 
reconnection. On the basis of a reconnection model, llsobe et al.\ (|2002[) showed 
that the energy release rate can be written as 

^ = SAJ r = ^B 2 c v m AJ r (1) 
at Ztt 

where, S is the Poynting flux into the reconnection region; B c , «i n , A r and 
/ r are coronal magnetic-field strength, inflow velocity, area of the reconnec- 
tion region and filling factor of reconnection inflow, respectively, llsobe et al.\ 
(|2002p have discussed the importance of filling factors in the estimation of en- 
ergy release rate, because not all magnetic-field lines reconnect. They suggested 
that it is best to use f r = 0.3. However, it is difficult to measure the coronal 
magnetic-field and a recourse to magnetic-field extrapolation is usually made. 
The reconnection area can be deduced from extreme ultraviolet (EUV) images. 
Estimation of inflow velocity (i>; n ) is another difficult issue, as there are very few 
direct observations for obtaining the inflow velocities (jYokoya ma et ai, 2001[ 



Narukage and Shibata, 2006). Following an indirect method to determine the 



inflow velocity from magnetic-flux conservation theorem, one can write 

B c Vi n = -BchroVribb = -Bphot^ribb (2) 

where, -Bphot and -B c hi are the photosphcric and chromospheric magnetic-field 
strengths, respectively, and u r ibb is the velocity of Ha ribbon separation. Here, 
we have assumed that the -B c hro ~ -Bphot- 

The inflow velocity Vi n can be thus calculated from Equation (|2|) using the 
values of B c , -Bphot and i; r ibb- But, accurate measurement of flare-ribbon separa- 
tion, f r ibb, is a difficult task. Different techniques have been used to measure the 
separation velocity of flare-ribbons, e.g., centre-of-mass motion, label matching, 
etc. Qu et ai\ (|2004bj) developed an efficient automatic technique based on 



component labeling and model matching, which is useful for flare forecasting. 
This method, however, provides only an average separation velocity of the flare- 
ribbons as a whole. To overcome this limitation, we have developed a simple 
technique which accurately measures the direction and velocity of each compo- 
nent of flare- ribbons. lAsai et aT1(|2006p have also studied this problem in detail 
for the event of 10 April 2001. They identified conjugate footpoints from cross- 
correlation by analysing the light curve at each Ha kernel, and tracked them to 
calculate separation velocity of flare- ribbons. They estimated the energy-release 
rate using magnetic-reconnection model ( |Asai et ai, 2004| from photospheric 
magnetic-field and Ha ribbon separation. They found that HXR bursts are 
formed around the peaks of energy-release rate. 
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3. Observational Data 

Wc have selected a large X17/4B flare observed in NOAA 10486 on 28 Octo- 
ber 2003 to demonstrate our technique. High spatial and temporal resolution 
Ha hltergrams obtained from the Udaipur Solar Observatory (USO), Udaipur 
(India) have been used for obtaining the chromospheric flare-ribbon separation 
velocity. These hltergrams were taken by a 15-cm aperture f/15 Spar telescope at 
a cadence of 30 seconds during the quiet phase, and 3 seconds in the flare mode. 
The spatial resolution of USO hltergrams is 0.4 arcsec pixel -1 . Photosphcric 
linc-of-sight magnctograms were obtained from GONG instrument having spa- 
tial and temporal resolution of 2.5 arecsec pixel -1 and one minute, respectively 
dHarvey et al, 19880 . 

Although we have used a rather large and complex flare in this study, the 
technique is useful for studying any two- ribbon flare event, in active regions 
or spotless regions associated with erupting filaments, which show separation 
with time, e.g., Ha flare-ribbons, Doppler, and magnetic features, or for that 
matter, flare-ribbons observed in other wavelengths including He I and He II. 
For example, |Maurya and Ambastha| (|2009p have studied the white-light flares 
of 28 and 29 October 2003 where flare-ribbons in multi- wavelengths, Doppler and 
magnetic features were analyzed. The technique can also be apply to events if we 
can extract the ribbons or features using mathematical morphological operations 
(MMO). Wc have listed some such events in Tabic [TJ 



Table 1. A list of two-ribbon flares studied using the technique 





Date 


Start 


End 


Class 


Position 


NOAA 












(lat-lon) 


No. 


1 


28 Nov 98 


04:54 


06:13 


X3.3/3N 


N17E32 


8395 


2 


17 Nov 99 


09:47 


10:02 


M7.4 


N18E17 




3 


17 Apr 02 


07:46 


09:57 


M2.6 


S15W42 


9906 


1 


28 Oct 03 


09:51 


11:10 


X17.2/4B 


S16E04 


10486 


5 


29 Oct 03 


20:37 


21:01 


X10/2B 


S15W02 


10486 


6 


20 Nov 03 


07:35 


07:53 


M9.6/2B 


N01W08 


10501 



4. Data-Processing Technique 

The raw Ha hltergrams obtained for our study require processing and correction. 
There are two main image processing steps to be carried out: i) pre-processing 
or data reduction and ii) post-processing or feature development to extract the 
feature of interest using MMO and other techniques described in the following. 

4.1. Pre-processing or Data Reduction 

We applied the following reduction procedures: From the large set of available 
images for the X17/4B flare of 28 October 2003, we selected the best images 
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taken in good "seeing" conditions. We applied (a) dark-subtraction to remove 
the signal due to dark current and (b) flat-fielding for equalization of the CCD 
pixels' response. 

Image alignment (or registration): Effects due to random tracking errors and 
image motion were removed by registration of the images. This is done in two 
steps: First, the images were manually registered using a compact sunspot, which 
was assumed as a fixed reference during the period of observation. Then, a 
second-step registration was carried out using a Fourier technique based on a 
cross-correlation method on the manually registered images to get co-aligned 
filtcrgrams registration within a sub-pixel accuracy. 

Intensity level normalization (or atmospheric correction): The average inten- 
sity level of solar images changes with time of the day, and also due to local 
effects of varying observing conditions caused by dust and clouds. To remove 
these effects, we carried out the following correction to obtain normalized images: 



where, I ln is the corrected input image obtained from previous steps, J m ; n and 
-fmax desired minimum and maximum intensity values in the output image I ou t, 
respectively. A corrected and normalized Ha filtcrgram is shown in Figure [2jt. 

We can set value of F min to and F ma >, to 255 for eight-bit images. But, the 
problem with this is that a single outlying pixel with either a high or low value 
can severely affect the value of F m j n or F max and this could lead to inappropriate 
scaling. Therefore a better approach is to first take a histogram of the image 
and then select F m - ln and F max at, say, the 5th and 95th percentile values in 
the histogram. That is, 5% of the pixel in the histogram will have values lower 
than F m i n , and 5% of the pixels will have values higher than F max . This helps 
in preventing the outliers pixels affecting the scaling significantly. 

4.2. Post-processing or Feature Development 

For our study, we are interested in flare-ribbons, which are usually the brightest 
features in Ha filtcrgrams. In order to select the regions of flare-ribbons in 
each filtcrgram, we set an appropriate threshold minimum intensity value. For 
example, a pixel value is set to zero, if it lies below 70% of the maximum value 
(c/., Figure (2b). The threshold was decided from a plage intensity. The flare 
intensity is taken to be greater than this threshold. The larger threshold will 
lose details of flare information while smaller value will include regions lying 
outside the flare- ribbons. In the case of lower (upper) threshold the flare-ribbon 
will be thicker (thinner) but will not affect the centroid of the ribbon from 
where the position measurement were carried out. Therefore, the calculation of 
flare-ribbon separation and velocity will not be sensitive to the threshold (see, 
Section [5]). It is, however, a difficult task to decide on the boundary of flare- 
ribbons for which one can use boundary-based methods. We have used a Sobel 
edge detector adopted in IDL™. Edges of flare- ribbons in Figure [2p, detected 
using this method, are shown in Figure [2]:. 
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Figure 2. Results of image processing steps applied to a typical Ho filtergram for the X17/4B 
flare event of 28 October 2003 at 11:03UT: (a) image obtained after the preprocessing, (b) 
region-growing, (c) edge-detected using Sobel method, (d) erosion of b, (e) dilation of b, (f) 
opening of b, (g) closing of b, (h) morphological closing of c, (i) small parts removed from h 
and (j) hole filled in i. Images b-j are obtained in negative for clarity. 



After obtaining the pre-processed images, we proceed to enhance and extract 
features of interest, e.g., flare-ribbons, using MMO, which is a powerful tool to 



extract the main feature of a digitized image (Gonzalez and Woods. 2008). In im- 
age processing, this tool was used for the first time bv ISerral(|1983j) to find the ge- 
ometrical structures in images. More recently, it has been used in solar physics for 
detection of flares ( Qu et a l. 2004a, 2004b|) , filaments ( Shih and Kowalski, 2 003; 
Qu et al., 2005| |Aboudarham et al., 2008|, su nspots (jZharkov et q/120051 I^OQgj) 
and prominences ( Fu, Shih, and Wang, 2007] ). Another technique based on neu- 
ral networks has also been used for the detection of solar flares by Fernandez 
Borda et al. ([2lM|) . 

The language of MMO is based on set theory. Sets in mathematical morphol- 
ogy represent objects in an image. A binary image is a complete morphological 
description of the image in 2D integer space (Z 2 ), where elements are either "0" 
or "1". The fundamental MMO are "erosion" and "dilation". If A and B are 
sets in Z 2 space then, "erosion" of A by B is defined as 



AeB = {z: {B) z C A} 



(3) 



In this study, A is the image (e.g., Figure 2b) to be eroded and B is the 3x3 
structuring clement. In erosion, smaller size of structuring element will eliminate 
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smaller components while larger size will eliminate larger components. Larger 
size will deform the shape of flare-ribbons which will change their centroid. This 
may affect the measurements of flare-ribbon separation and velocity. A more 
description about the structuring element with different MMO can be found in 
IGonzalez and Woodsl ([2008]) . 

The "erosion" shrinks the components of an image. It consists of replacing 
each pixel of an image by the minimum of its neighbours (Figure EH). The 
"dilation" of A by B in Z 2 space is defined as 

A®B = {z: {B) z nA^4>} (4) 

The "dilation" expands the components of an image. It consists of replacing each 
pixel of an image by the maximum of its neighbours (Figure [2J3) . There are two 
other morphological operators: "opening" and "closing". The "opening" gener- 
ally smooths the contour of an object, breaks narrow isthmuses, and eliminates 
thin protrusions. The "closing" also tends to smooth sections of contours but, 
as opposed to "opening", it generally fuses narrow breaks and long thin gulfs, 
eliminates small holes, and fills gaps in the contour. The "opening" of a set A 
by structuring element B is defined as 

AoB = {AQB)®B (5) 

i.e., opening of A by B is the erosion of A by B, followed by a dilation of the 
result by B, Figure Hf. Similarly, "closing" of set A by structuring element B is 
defined as 

A*B=(A®B)QB (6) 

i.e., closing of A by B is the dilation of A by B, followed by an erosion of the 
result by B (Figure Hk) . We applied morphological "closing" to the image d2t) 
to erase the gaps and smooth the contours. A binary image after this operation 
is shown in Figure [2ji. The images after applying MMO consist of many small 
features. These are removed using "region labeling" methods. The region labeling 
method gives labels to each components of an image. We count the number of 
pixels in each label and remove those component by "0" , which have values less 
than a set threshold (Figure [21) . 

After removing small components from images obtained using the above 
method, we find small holes, or background regions surrounded by a connected 
border of foreground pixels within selected components, i.e., large area compo- 
nents as seen in Figure^. In fact, the holes in between the flare-ribbons are part 
of the flare-ribbons, created due to edge detection by the Sobel operator. Since we 
are interested only in the edge, area, etc. of flare-ribbons, these holes should be 
removed by "1" from binary images. For this purpose, we have used an automatic 



procedure based on morphological reconstruction (Gonzalez and Woods, 2008). 
This procedure uses "geodesic dilation" . This is described as follows. 

Let I(x, y) be a binary image. We form a marker image F, which is "0" 
everywhere except at the image border, where it is set to 1 — /; i.e., 

„, s _ J 1 — I(x, y) if (x, y) is on the border of I ,_, 
*{x,y) - \ o otherwise [ ' } 
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Then, the binary image with holes filled will be, 

H = [R%(F)] C . 



(8) 



Superscript c represents the complement and [i?^(i 7 ')] c is "geodesic dilation" of 
size n of the marker image F with respect to the mask image I c . The "geodesic 
dilation" is defined by. 



(9) 



where, R^) (F) = (F © B) n I c . n represents the "intersection" operator and B 
is the structuring element as described above. Image after removing the holes is 
shown in Figure [2- 



5. Determination of the Flare-ribbon Separation 

Once the flare-ribbons are extracted, we can estimate their separation from a 
refrence, such as the magnetic neutral line. Some questions that arise in the 
determination of flare-ribbon separation are: Can we use a straight line repre- 
senting the neutral line passing between the flare-ribbons? Do the flare-ribbons 
move perpendicularly to the neutral line? Do the flare-ribbons on either side 
move with the same velocity? From the visual inspection of the Ha flare-ribbons 
(c/., Figure EK), we find that the ribbons are nearly parallel to each other, but 
they are curved in shape. Therefore, we can not draw a straight line representing 
the neutral line. Also, various components of the flare-ribbons are curved over 
different directions. Therefore, the direction of neutral line is "important" to get 
the direction of motion of a given part of the flare-ribbon. 

Figure [3] shows a cartoon of two ribbon flare to illustrate the method used 
here for distance measurement between flare-ribbons and the neutral line (NL). 
Ri and R2 represent two ribbons of the flare, and the dotted line NL represents 
the magnetic neutral line. Solid lines, drawn perpendicular to the neutral line, 
mark the directions of motion of different portions of ribbons, assuming that the 
apparent motion of flare-ribbons is perpendicular to the neutral line. Therefore, 
the velocity derived for a given point on the flare-ribbon may be the same in 
magnitude, but it may not necessarily be so in direction. It is important to con- 
sider this issue appropriately in the analysis. Further, it is obscrvationally known 
that the velocity of flare-ribbons is affected by the strength of the magnetic-field 
in the region of flare-ribbon motion. However, it should be noted that this is 
only an apparent motion which may not represent the real motion of the flare 
arcade. 

It is well-known that the overlying arcade has shear, which decreases with 
increasing distance from the magnetic-polarity reversal boundary and with in- 
creasing coronal height. Thus it follows that footpoints of the reconnected flare 
loops will move not only from the polarity reversal boundary, but also along it. To 
address this issue, we need to obtain 3D topology of the magnetic-reconncction 
region. We should also know the conjugate points joining the flare-ribbons, 
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which is not easy. Some information about this can be obtained from magnetic- 
field extrapolation and post-flare loops seen in EUV images. For example, to 
find the conjugate points in Ha flare-ribbons, lAsai et a l. (2006) divided the 
regions of flare-ribbons lying over magnetic polarities opposite to each other 
into fine meshes. Then, they drew light curves of total intensity for each box 
in both meshes, and identified the highly correlated pairs (conjugated pairs) by 
calculating the cross-correlation functions. This is a good approach to find the 
conjugate points, but for a complex region, field lines may be highly sheared 
in small scales which poses difficulty in identifying these points. The magnetic- 
field extrapolation may be a more accurate approach to get much closer to 
the conjugate pairs. But, again, for complex regions errors in magnetic-field 
extrapolation pose difficulty to identify the exact conjugate points. Therefore, 
we have taken a simple assumption of flare-ribbon motion perpendicularly to 
the neutral line. 

In view of the above-mentioned issues, we have implemented an algorithm 
in IDL™ for detection of flare-ribbon separation based on an automatic tech- 
nique. It calculates the ribbon separation measured from the neutral line at 
specified points of the neutral line. The algorithm assumes that the ribbons 
move perpendicularly to the neutral line. The main steps of the technique are: 

• Neutral line passing between the flare-ribbons: We overlaid contours of co- 
spatial photospheric line-of-sight (LOS) magnetograms over the correspond- 
ing Ha filtergrams observed around the time of the flare. The neutral line 
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(NL) drawn in green colour is obtained by fitting an appropriate polynomial 
over several points within the region of magnetic-field lines drawn at ±5 
Gauss levels (Figure HI left panel). 




Figure 4. Left panel: Ha filtergram of NOAA 10486 overlaid with the contours of longitudinal 
component of magnetic-field showing the X17/4B flare of 28 October 2003 at 11:03 UT. Solid 
and dashed contours represent positive and negative polarities, respectively. Dotted curves 
represent the polarity reversal lines. A part of neutral line NL passing through the flare-ribbons 
is highlighted in green colour. Right panel: Extrapolated magnetic-field lines plotted over the 
GONG LOS magnetogram using the IDL package MAGPACK2. 

• Direction of motion of the flare-ribbons: For simplicity, we assume that 
flare-ribbons move perpendicularly to the neutral line. Let NL be the neu- 
tral line (c/., Figure [3]) passing through the flare- ribbons R\ and R2, and 
P be a point on it. AB is tangent at P on NL, and CD is perpendicular 
to AB at P. Let Q be the centroid of the part of flare-ribbon R2 which 
lie along PD and follows path PD. It is required to find the direction of 
PD with respect to the horizontal or x-axis. For this, we find the gradient 
(say, m p ) at point P derived by Lagrange interpolation method. Therefore, 
the gradient of line PD will be m' p = —\/m p . A similar procedure can be 
carried out to find the direction of motion of a flare-ribbon at any desired 
point along NL. To obtain the coordinates of points on a line, say PD, we 
have used a simple algorithm as described in the Appendix. 

• Distance d t of a point Q(xk,yk) on the flare-ribbon at time t from point 
P(x p ,y p ) on the neutral line is given by 



dt = yj {xk ~ Xp) 2 + (Vk - Vp) 2 
and the direction of motion is 

8 p = ta?i _1 (mp) 

Using the calculated flare-ribbon distances from the neutral line, one can find 
the velocity simply by taking their time derivatives. But this would be noisy due 
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to the errors in the measured distances using above procedures. Therefore, an 
appropriate function is used to fit the measured distances so that the errors are 
reduced. 



6. Ha Flare-ribbon Separation and Energy-release 

The X17/4B event of 28 October 2003 was a two-ribbon white light flare that 
occurred in the super active region NOAA 10486. An Ha filtergram for this 
event at 11:03 UT overlaid with contours of longitudinal magnetic-field shown 
in Figure [4^. Solid and dashed lines represent positive and negative magnetic 
polarities, respectively, and the dotted lines represent magnetic neutral lines. 
The neutral line (NL) passing through the flare-ribbons is highlighted in green. 
Straight lines drawn perpendicular to the neutral line mark the directions of 
motion of different parts of flare-ribbons from the corresponding points Pi (i = 
1, . . . , 6) of the neutral line. Here, we have selected six points Pi over the neutral 
line and corresponding directions of flare-ribbon motion are shown. The center 
of the two edges of a flare-ribbon is followed along the line passing through P 
to measure the flare-ribbon separation from corresponding points (i.e., p) on 
NL. The measured distances corresponding to different components are shown 
by "x" symbol in the panel (a) of Figures ©-([7]). 

Velocities corresponding to six selected parts of the flare-ribbon were obtained 
by taking time derivatives of the fitted distances. We found Boltzmann-Sigmoid 
to be the best fit for this case ( jMaurya and Ambastha, 2009 ). The fitted dis- 



tances are shown by solid lines passing through the measured points in the 
panel (a) of Figures {5]) — (0). Corresponding velocities are shown by dotted lines 
in the same panels. 

6.1. Magnetic-reconnection and Energy-release 

Using the derived separation velocities for the flare-ribbons, we can estimate 
magnetic-reconnection rate and energy released during the solar flare. Magnetic- 
reconnection theories predict how fast reconnection can occur. The speed of 
magnetic-reconnection can be expressed in terms of the inflow velocity v- m or 
the dimensionless ratio of Vi n to the Alfven velocity (va) in the inflow region. 

V\r> 

Mi = — 

where, the Alfven velocity is given by va = B c / '^Anp and p is density near the 
reconnection region, which can be taken as a free parameter from an atmospheric 
model. Therefore, the reconnection rate using Equation {2]) can be written as, 

Mt = ^BphotV^ (1Q) 

An alternative measure of the reconnection rate is the electric-field strength 
in the RCS. It shows how violently the magnetic-reconnection progresses. It is 
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Figure 5. Top and bottom panels correspond to the measurement of different parameters 
along the lines P± and P2, respectively (c/., FigureQ. (a) "x" represents the measured dis- 
tances, and the solid curves passing through these are the Boltzmann-Sigmoid fitted distance 
profiles with time. Dotted lines represent the velocity v(t) profiles derived by taking time 
derivative of the fitted distances, (b) Magnetic-flux B p ^ ot at the points used for distance 
measurement, (c) Reconnection rate (<J>) (solid line) and Poynting flux (S) (dotted lines). 
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Figure 6. Top and bottom panels correspond to the measurement of different parameters 
along the lines P3 and P4, respectively (c/., Figure |4j. (a) "x" represents the measured dis- 
tances, and the solid curves passing through these are the Boltzmann-Sigmoid fitted distance 
profiles with time. Dotted lines represent the velocity v(t) profiles derived by taking time 
derivative of the fitted distances, (b) Magnetic-flux B p hot a t the points used for distance 
measurement, (c) Reconnection rate ($) (solid line) and Poynting flux (S) (dotted lines). 
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Figure 7. Top and bottom panels correspond to the measurement of different parameters 
along the lines P5 and Pg, respectively (c/., Figure [4} . (a) "x" represents the measured dis- 
tances, and the solid curves passing through these are the Boltzmann-Sigmoid fitted distance 
profiles with time. Dotted lines represent the velocity v(t) profiles derived by taking time 
derivative of the fitted distances, (b) Magnetic-flux B p hot a t the points used for distance 
measurement, (c) Reconncction rate ($) (solid line) and Poynting flux (S) (dotted lines). 
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defined as the reconnected magnetic-flux per unit time and is expressed as, 

$ = B c v in 

From Equation it can be written as 

$ = £ p hot«ribb ( n ) 

The rate of magnetic energy-release, i.e., energy released per unit time during a 
solar flare, is the product of Poynting flux and the area of RCS, given by Equation 
([1]) . If we assume that the area of magnetic- reconnection is fixed during the flare 
(as assumed bv lAsai et al.\ l2004|) . the energy-release rate will be proportional 
to the Poynting flux only and hence independent of the filling factor / r . Using 
Equations fl} and (|2|), Poynting flux can be written as, 

S = — -Bc-Bphot^ribb (12) 
Z7T 

The amount of energy reaching at a footpoint is not the total energy because 
a part of the reconnection energy will go outward and only the rest will come 
downward to the solar surface. Further, the energy coming toward the solar sur- 
face will be distributed into two footpoints. Therefore, the energy at a footpoint 
can be written, ( |Asai et at, 2004[ ) as, 

S = — -BcSphot^ribb (13) 
Z7T 

where, factor e is (0 < e < 0.5). In the above equation, the coronal magnetic-field 
(B c ) is yet to be determined. The numerical calculation of Mi and S requires 
B c . To replace B c by the photospheric magnetic- field -B p hot, lAsai et ai\ (|2004|) 
considered that B c = a_B p hotj where a is a constant to be measured. The Poynting 
flux then can be written as, 

S = ^Bphot^ibb. (14) 

From the magnetic-field extrapolation (using IDL™ package MAGPACK5£], 
Sakurai| 1982)) lAsai et all (|2004|) found the ratio of coronal to photospheric 
magnetic-fields (i.e., a) m 0.2 for the X2.3 flare event of 10 April 2001. Using 
this procedure, we extrapolated the GONG LOS magnctogram (c/., Figure |H 
right panel) to find the constant a for the X17/4B flare event of 28 October 2003 
at 11:10 UT. This turns out to be « 0.09, which is an average value over the 
whole region, and is smaller in strong field area. 

To calculate the magnetic-reconnection [Equation pip ] and energy-release 
rate [Equation (fT4|) ]. we require the strength of the photospheric magnetic-field. 
For this, we estimate the photospheric magnetic-field (Bphot) from GONG LOS 
magnctogram at the center of the two edges of a Ha flare-ribbon corresponding 



http: / /solarwww. mtk.nao.ac.jp/sakurai/en/magpack2.shtm I 
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to i> r ibb measurements. The estimated -Bphot is shown in Figure ((5]b) - (JTJd) . The 
reconnection parameter ($) and energy-release rate (S) are shown in Figure 
(03 ) - dZt) with solid and dashed lines, respectively. 

|Xie, Zhang, and Wang| (|2009p studied flare-separation speeds and reconnec- 
tion rate of the 10 April 2001 solar flare and reported a weak negative correlation 
between the ribbon-separation speed and the longitudinal magnetic-flux density 
as found here. However, contrary to the event studied by them, we do not find 
similarity in time profiles of the magnetic-reconnection rate and the flare-ribbon 
separation speed during the evolution of the X17/4B flare. In fact, we find 
that the magnetic-reconnection rate is better correlated with photospheric LOS 
magnetic-field strength. 

The average Poynting flux over the entire RCS during the flare is found to 
be of the order of 10 9 . This flare lasted for over 10 4 seconds. The approximate 
area of RCS estimated from EUV is 10 19 cm 2 . Therefore the energy-release is 
of the order of w 10 32 ergs, which is sufficient to produce the HXR sources (c/., 
Figure [TH). 

7. Summary and Conclusions 

We present an automated technique to extract the flare-ribbons on the basis 
of mathematical morphology. Using the assumptions that flare-ribbons separate 
perpendicularly to the neutral line, we developed an automatic algorithm to com- 
pute the motion of flare-ribbons. Using flare- ribbon separation and photospheric 
magnetic-field, we computed the magnetic-reconnection rate and energy-release 
rate during the large solar flare X17/4B in NOAA 10486 observed in Ha at the 
USO. 

This event occurred in a complex active region showing both extended ribbon 
structures as well as several fragmentary kernels. The flare was associated with 
filament eruption, two-ribbon separation, and a very fast CME. Temporal evolu- 
tion of separation between different portions of the two main ribbons is derived. 
Nearly all of these show a rapid increase of velocity reaching a maximum in the 
range 20-60 km s _1 in the initial stage, i.e., 10:58-11:05 UT. Thereafter, during 
the next ten minutes, the separation speed decreased at different rates depending 
on the magnetic-field structure and strength in the region. Interestingly, the 
magnetic- flux showed a rapid decrease/increase in some locations at the time of 
maximum velocity of ribbon separation. 

Magnetic-reconnection in corona, as represented by the electric fields (E c ) 
in the reconnecting current sheet, has been estimated from the measured rib- 
bon separation speed and magnetic-fields obtained from GONG. Electric fields 
reaching up to 35 kV m^ 1 were found in some locations. The evolution of these 
parameters provides evidence that the impulsive flare-energy release is indeed 
governed by the fast magnetic-reconnection in the corona. The magnetic-energy 
release rate (i.e., Poynting flux) has also been deduced for various parts of the 
flare. An average Poynting flux during the flare is found to be of the order of 
10 9 , which gives the energy-release of the order of f» 10 32 ergs over the entire 
flare duration of over 10 4 seconds. 
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In summary, we have found that: i) the entire flare-ribbon does not move 
with the same velocity, but different parts of flare-ribbons move with widely 
varying speeds in different directions, ii) the flare-ribbon separation is observed 
to be decelerated in the regions of strong magnetic-fields, and Hi) magnetic-field 
reconnection (i.e., electric-field strength) and energy-release rate (i.e., Poynting 
flux) increase with increasing separation velocity. 

A complete and accurate flare-energy release calculation is beyond the capa- 
bility of the available observational and theoretical techniques. The 3D topology 
of active regions is only now beginning to be observed by recent projects such 
as the Solar TErrestrial RElations Observatory (STEREO) spacecraft. The aim 
of the present study is to provide a relatively simple, automated approach to 
estimate the approximate level of energy-release in two-ribbon flares. 
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Appendix 

Coordinates of the points of a straight line: Suppose, we have to find coordinates 
of the points of a straight line RS, with points R and S having coordinates 
(xi,yi) and (2:2,2/2), respectively. The distances between abscissa and ordinates 
of the two points will be, n x — abs(a;2 — %i) and ny = abs(j/ 2 — Hi), respectively. 
Let, n = max(n x , n y ). Now 

if nx > ny then 

if xl = x2 then xi = xl and yi = yi+ i 
if xl 7^ x2 then yi = m[xi — xl) + yl 

else 

xl = x2 and yl = y2 

where, m = (1/2 — yi)/(x2 — x\) is the gradient of line RS and i = 1, n are the 
indices of the coordinates of points of the line RS. 
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